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Abstract 

This paper discusses the elastic behavior of a single polyelectrolyte chain. 
A simple scaling analysis as in self avoiding walk chains are not possible, be- 
cause three interplaying relevant length scales are involved, i.e., the Debye 
screening length and the Pincus blob size. Therefore a self consistent compu- 
tation of an effective variational propagator is employed. It is shown that the 
elastic force as a function / of the distance R behaves as f on R for small /. 
For larger forces we find a new regime, characterized by deformations larger 
than a computed electrostatic "blob size". These results are supported by 
simulations and intuitive physical arguments. 

PACS: 05.20.-y, 36.20.-r, 61.41.+e 

short title: Elasticity of polyelectrolyte chains 
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1 Introduction 

In contrast to neutral polymers polyelectrolytes are bearing ionizable groups, that 
are able to dissociate into charged polymer chains and small counter-ions. The 
interest in polyelectrolytes reaches back to the early days of polymer science (see 
e.g. [0]), due to their fundamental importance in biology, biochemistry as well as 
in industrial applications. Proteins, nucleic acids and superabsorber materials are 
only a few examples out of the wide range of practical interest in polyelectrolytes. 
Nevertheless they belong to the least understood systems in macromolecular science 
[0]. The reason why they are much less understood than for example neutral poly- 
mers lies in the difficulty to apply renormalization group theories and scaling ideas 
to systems in which long range (i.e. coulomb) forces are present, which means that 
the range of the forces is of the order of the polymer size. 

For many applications knowledge about the elastic behavior of polyelectrolytes 
is of fundamental importance. One of them is, for example, to understand the 
elastic behavior of superabsorber materials, i.e., polyelectrolyte networks and their 
thermodynamic properties. Most important amongst them is naturally the swelling 
behavior. 

The elastic and the thermodynamic properties of classical networks formed by 
crosslinking of neutral polymer chains have been considered by simple and successful 
theories |^ . The success of these theories lies in the fact that most of them are based 
on the assumption that the elasticity of the entire network is roughly given by the 
elasticity of a representative single chain in the cross-linked network (see for example 
[^], for reviews). Many theories including those which consider entanglements. 
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work very well for dry networks, i.e., networks that are not swollen by a good solvent. 
The case where the networks are at equilibrium swelling at the so-called c* point, 
when the blob size is precisely identical to the radius of gyration of the mesh |^ is 
especially interesting. In this case the network elasticity can be computed by the 
knowledge of the conformational behavior of a single swollen chain 

Naturally these simple theories are suitable for many purposes but are, of course 
on the other hand, unable to answer difficult questions such as the liquid solid 
transition, etc ||, |^, but this is not the subject of the present paper. We would 
like to address another important issue of the subject of network elasticity. In all 
theories the effects of interactions have been neglected. In the classical theories it 
was assumed by purpose of mathematical simplicity that any interactions between 
the chains to do not contribute to the elastic properties. This is most likely valid 
for neutral polymer networks, because in dense systems the excluded volume forces 
between the chains are screened to a large extend. In polyelectrolyte networks this 
is definitely not the case, because strong interactions between monomers rule the 
conformational behavior . Indeed little is known about the interplay of elasticity 
(conformation) and strength of the interactions in the theory of network elasticity. 

To get a first insight and physical feeling for the interplay of these effects, it is 
useful and legal to consider single chain elasticity in weakly charged and fiexible 
polyelectrolyte chains, although the simple generalizations as carried out for neutral 
networks are not promising. On the other hand simple scaling considerations as 



done for excluded volume chains [y, |Tl| cannot be employed a priory, because several 
independent length scales are involved. The first is (apart from the chain length) the 
range of the electrostatic interaction. For simplicity we use here a screened Debye 



potential V{r) oc l/rexp(— Kr), where Ad = 1/k defines the screening length of the 
electrostatic interactions. The other scale, in neutral SAW polymers introduced as 
Pincus blob size via the relation / = k^T/^p, defines the elastic properties as a 
function of the excluded volume interaction by the assumption that within scales 
the chains shows excluded volume behavior. / is the external force acting at the 
chain ends. (Since we assume / = constant, the same force acts on each segment 
along the chain). Similar treatment is here not possible, because the interplay 
between Ad and ^p is unknown. It has been shown only recently, how a field theory 
can be set up for the critical behavior, i.e., chain length ^ oo of undeformed 



polyelectrolyte chains [|T2 



As a first step we calculate the elastic response of a single polyelectrolyte chain 
in solvent. Therefore we apply an constant external force on the ends of the chain, 
which corresponds to the application of an external field. Here we follow the simple 
idea in classical rubber theory, but determine the new propagator of the extended 
interacting chain. To do so, we employ a variational principle and calculate the 
effective propagator of the chain, which allows statements about the influence of 
conformation and interactions on the elasticity. Although this procedure seems 
to be oversimplified for several reasons. The first is that it cannot be expected 
that polyelectrolyte networks can be described by effective single chain models as 
in networks with short range interactions. The second reason is that fluctuations 
are expected to play an important role. Thirdly, the Debye-Hiickel approximation 
for the electrostatic potential does not seem to be sufficient. This can be seen in 



13] where the Debye-Hiickel potential simply appears by a random phase approx- 



imation in a Gaussian density functional theory. Moreover there is great evidence 



from numerical simulation that the Debye-Hiickel approximations fails under certain 



circumstances . However, the use of the Debye approximation together with the 
assumption that the non interacting chain is Gaussian limits the model to weakly 
charged and flexible polyelectrolytes. Nevertheless, despite these drawbacks, we 
expect first principal result from the theory below. As such a result we derive the 
elastic modulus of the single polyelectrolyte chain and compare it the the well known 
cases, i.e., free chain and excluded volume chain. As a side remark we mention that 
the same method can be used to re-derive the Pincus blob model for pure excluded 
volume chains. In this case the renormalized propagator can be determined in a 



field theoretic A — d - expansion fTS 



The starting point of the calculation is the Green function of the free chain 
without interactions. The Green function can be viewed as the Fourier - Laplace 
transform of the distribution function for the chain ends. It will be shown that 
the force in the system is treated through the analytic continuation of the Fourier 
transformed Greens function to the complex plane as will be shown in section two. 
After having introduced a field theory the problem is mapped on a Gaussian field 
theory with a propagator that formally in the Fourier space can be written down 
exactly by making use of the proper self energy. According to the well known 
Feynman variational inequality the sum of the Gaussian free energy and the mean- 
value of the interacting potential has to be minimized with respect to the proper 
self energy, which is our variational parameter. This leads to a non-linear integral 
equation for the proper self energy, which can be solved systematic approximations. 

The result of this consideration is that the chain is stretched in the long ranged 
limit proportional to the applied force parallel to the force and is pushed proportional 



to the force perpendicular to it. This result is in good agreement with the simulation, 
that is also presented in this paper. Moreover we support the computations by 
intuitive physical arguments. The variational calculation and simulation indeed 
support a scaling picture of the deformation of the polyelectrolyte chain. 

The paper is organized as follows. In the next section we present the mathemat- 
ical model and provide the definitions. In section 3 we introduce a field theory and 
calculate the variational equation for the proper self energy. In section 4 this varia- 
tional equation is solved approximately. The results of the analytic calculation are 
presented in section 5. The results of the computational simulation are presented in 
section 6. Section 7 is provided with some scaling considerations. 



2 Model 

Let us first introduce the standard model which is employed here. Since we restrict 
ourselves to flexible chains which are weakly charged the Edwards model is the ap- 
propriate tool. The starting point is thus the continuum version of the dimensionless 



Edwards Hamiltonian 16 



3 r^o , fdvV . , „dr 



+ — / ds ds' ^ \ ' ^ \ ^'^ , (1) 



2 Jo Jo \r{s) — r(s') 

where r(s) represents the chain conformation in three dimensions as a function of the 
contour variable s, b = /ATreoerksT is the Bjerrum-length, (3 is {ksT)^^, where ks 
is the Boltzmann constant and T denotes the absolute temperature. I is the Kuhn 
segment length, z is the monomer charge in units of e, eg is the dielectric constant 
and er the relative dielectric constant. Nq stands for the bare number of monomers 



on the chain, f is the external force and k ^ denotes the Debye-Hiickel screening 



length. 



The correlation function can be calculated in terms of a path integral |T6|, |T^, ^ 
as follows: 



G (r, iVo; f) = / Vv{s) exp {-(3HE[r; f]} 

Jr(0)=O 



(2) 



Its Fourier transform is defined by 



G(k,A/'o;f) = j d^r exp{-ikr}G (r, iVo; f) 



(3) 



The averages of the force-size relationship (Ry) and (Rj_), where Ry denotes the 
parallel component with respect to f and R^ is the corresponding perpendicular 



part, are then readily calculated by the general formulae 



(R^)[f] 



d^/dkl G(k, iVo; f ) 



G(k,iVo;f) 



(4) 



k=0 



and 



(Rl)[f] 



ELi9V5A;2G(k,iVo;f) 



G(k,iVo;f) 



(5) 



k=0 



By analytic continuation of the Fourier space to the complex plane, the correlation 
function (^(k, A'^q; f) can also be written as the zero-force correlation function G{k — 
ipf, No;f = 0). Substitution of Eq. |^ into Eq. || yields: 

r „ fr{No)=r 

G(k,No;i)= dh exp{-ikr} / Vr{s) exp {-/Jf/^ir; f]} (6) 

J Jr(0)=O 

For constant force f Eq. (^) can be rewritten as: 



G(k,iVo;f) 



l'r{No)=r 

dh exp{-i(k - z/3f)r} / Vr(s) exp {-PHeIt; f = 0]} 

Jr(0)=O 



^(k - ipf, No; f = 0) 



(7) 



Consequently, to get results for (R|) and (Ri), we only have to calculate G{'k,No) 
and continue the first argument of G to the complex plane. 



3 Field-theoretical formulation 

The Laplace transform of G{k, No) with respect to A^o is defined by 

G'(k,/xo)=/ diVoexp{-/ioiVo}G(k,iVo) (8) 

The function G(k,fio) can be calculated by the introduction of de Gennes' zero- 
component field theory (see for example [PD|]) 

G{k,i,o) = lim / D^^i(k)Vi(-k)exp{-/5i7[V^]} (9) 

Here the field theoretical Hamiltonian Hltp] is given by 



o 



m) (10) 



+ ^ I V^(ki)7/;(k2)[/(ki + k2)(5(ki + k2 + ks + k4)^(k3)V^(k4) 

O Jki,k2,k3,k4 

where /j^. is an abbreviation for / (i^k/(27r)^ and U denotes the Debye-Hiickel poten- 
tial in units of (3~^. In the Fourier space (^(k, yUo) can be written exactly as 

G'(k,/io)= + + S(k)'| (11) 



where S(k) denotes the proper self energy. We now consider an approximate corre- 
lation function ^(k, /xq) with an approximate proper self-energy A^(k). 

^(k,/io)= (^/io + ^A;2 + ^(k)^ (12) 

Defining the Hamiltonian Ti by 

= \ Jj{-k)g-\k,f,o)m) (13) 



^(k, /io) can be calculated in the following way 

^(k,//o) = lim / V^M^)M-^)eM-P^m (14) 
In this notation the well-known Feynman inequality is given by: 

F <T+{H-n)H (15) 

where 

^ /W...exp{-^H} 

"-0 /I)^ exp{-/37^} 

is the mean-value and T the free energy with respect to Ti. The right hand side of 
the inequality (|15D has to be minimized with respect to A^. ^ and (if — 'hOj-n can 
be written in terms of the correlation function ^(k, yUo): 

= - In Z = -^y ^ In ^(k, /io) (17) 

As can be shown easily the second term of the right hand side of inequality ( |T5| ) is 

/-, 2 2 2 
U 2 [ g(ki,/io)g(k2,/io) 

+ Txhz nV / — ^— 18 

The general minimization condition reads 

(^+(/J-7^)«) = (19) 



5M(q) 

where 5/(57W(q) denotes the functional derivative with respect to 7W(q). After 
inserting Eqs. ([17|) and (|18]) into Eq. (|19]) one obtains 

2'Kbz'^n 



^^"^^ A^o + |P + A<(k) 



+ 47r6z2 /■ J ^ ^ (20) 

7k («:2 + (q + k)2) (^0 + f + -M(k)^ 
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This is a non-linear integral equation for Ai{q), which in the following has to be 
solved approximately, since the exact solution is unknown. At this point it should 
be stressed, that Eq. ( ^0]) represents the well known Hartree approximation. 

Another important and useful point is, that the exact proper self-energy S(k, /iq) 
is less than or equal to the approximate proper self-energy J\4(k, (Iq). This can be 
shown as follows: First of all we introduce an infinitesimal auxiliary real field h{r) 
in the field theoretical Hamiltonians H (Eq. (0)) and H (Eq. (|I^). 

PH[^f] ^ pHliJ] + I /i(k)^i(-k) (21) 

An analogous extension has to be done for 7Y, where h{\s) is the Fourier transform 
of hijc). Using this Hamiltonian the exact free energy F becomes a functional of the 
auxiliary field h. Thus: 

F[h] = - In (/ PV^e-'^^-/.^^«^^(-'^)) (22) 

Evaluating Eq. (|2^) for small h yields 

F[h] = F[0] |/.(k)|^G(k,/io) + 0{h') (23) 

The validity of Feynman inequality is unaffected by the introduction of the auxiliary 
field h and becomes after having neglected 0{h'^) 

F[0] - i 1^ |/i(k) 1^ G(k, ^o) < -^[0] + {H-'H)n-\ \Hk)\' ^(k, /xo) (24) 

Inequality (^3]) can be rewritten as 

i \h{k)f (^(k, /.o) - G(k, /.o)) < m + {H- H)n - F[0] (25) 

In the limit n ^ the right hand side of Inequality (|25|) vanishes. Due to the fact 
that inequality (|^) holds for any field h, we get 

G(k,/io) >^(k,/io) (26) 
11 



which is equivalent to 



M{k) > S(k) 



(27) 



Thus A^(k) is proven to be an upper bound for S(k). 



4 Approximate solution for the proper self- 



energy 



To do the exphcit calculation let A^r(q) be A^(q) — A^(0). Then A4r{<i) is given 

by 



1 



1 



/i + fP + 7W,(k) 



K2 + (q+k)2 /s:2 + k2_ 
where = fiQ + A^(0). In order to simplify the integral in Eq. 
following approximation, which is valid for small k: 



:28) 



we make the 



MM) 



|k|>K 



(q+k)2 



1 

k2 



fx + ^}k^ + Mr{k) 



(29) 



Equation ( p9D could be solved by means of an iteration procedure following the 



scheme 



q) 



|k|>K 



(q+k)2 



1 

k2 



J 



(30) 



q) 



|k|>K 



Airbz'^ 



q + k)2 k2j ^+ilk^ + Mr'\k) 



(31) 



but we show below that in the variational technique this procedure is not necessary. 
En effet, it is show right below, that the one loop renormalization agrees with the 



first order perturbation in the limits we investigate. 
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It can be seen from Eqs. (^) and ( ^iD that Air{^) is actually a function of the 
dimensionless parameters q/ and kI. Denoting Q = ql/nl and K = kl/nl Eg . (|28D 
becomes 



"'|K|>1 



1 1 



(Q + K)2 K2 



(32) 



In the limit of small kI the validity of the ansatz 

MriKd) = aK^K^f + 0{{KkIY) (33) 



can be checked using Eqs. (^0]) and (31). Consequently, the only remaining task is 
to calculate a from Eq. (|3^ ) selfconsistently. Therefore we introduce in the integral 
of Eq. (^) spherical coordinates [K^ i), cp) and perform the integration over cp and 
1?. The second derivative with respect to Q at Q = yields the following equation 
for the coefficient a: 

We perform the integral by only taking into account the most singular term with 
respect to k. The result of this calculation is according to Eq. ( p3D 



a 



(35) 



It is important to note that Eq. ( pS]) coincides exactly with the first-order term in 
the perturbation expansion. Therefore the higher order terms within the Hartree- 
approximation do not contribute to the coefficient a. This underlines the quality of 
the first-order approximation. 

The constant term of the approximate proper self-energy A^(0) is given by the 
expression 



TT 



= / — — 36 
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Since /io is a finite number greater than zero, the main contribution to the integral 
in Eq. (|36|) comes from large k. Eq. ( ^9]) shows, that limg^oo -^r(q) = —-^(0), 
because of the fact, that the first term in the brackets of the integrand can be 
neglected in the considered limit. Then the right hand side of Eq. (pOf ) is exactly 
minus the right hand side of Eq. (pG]). Thus for large k the integrand does not 
depend on A4 and therefore on k. Consequently, /i = fio + ^A{0) contains no 
singularity for small k. 

Note that we have confined ourselves to small values of q so far, i.e., such values 
for which q <^ n is satisfied. This makes sense, if we only consider the end-to-end 
distance of the chain without external force. In this paper, however, we introduce 
an external force on the chain. Therefore the restriction g -C k may be to strong 
according to Eqs. @ und (|^). 

Starting from Eq. (|30|) we obtain 



k|>K 



(37) 



The integral in Eq. ( ^7\ j can be calculated in a double expansion. The result is 



-M«(q) = 



I / q \ 2n+2 



2 (2n + l)(2n + 3) [J 
If one neg lects O ({KH'^/^f^) Eq. (H) becomes 



A<(l)(q) = 5!!^ifi_!i),nfi±4!iU^^i (39) 

^ ' Atx q) \l-q/K) 2n jj, ^ ^ 

As can be seen from Fig. (1) the quadratic approximation of Ai^^^ works well even 
for values such that q ^ k. It is very important to notice that Ai^^^ has an imaginary 
part as soon as g > k, which gives also the limit of the force /, since the the vector q 
contains the external force as third component. Indeed we need here the condition 

Pf/K < 1. 
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Figure 1: Graphical illustration of the validity or our quadratic approximation for 
the proper self energy: quadratic approximation (dashed function), exact function 
according to Eq.(38) (bold-faced line) 



An inspection of Eq. (37) shows that the appearance of the imaginary part is 
clearly an artifact. Therefore terms of 0{nl/ ^/Ji) cannot be neglected even if nl/ ^JJl 
becomes very small. For q > k, the quantity nl/^/Ji cannot be used as a small 
parameter within the perturbation expansion. As a consequence we expect a new 
regime which will be detected also in the simulation and the scaling theory below. 

5 Results 

Inserting the approximate result for the proper self-energy S (k) into Eq. (|TI|) yields 
an explicit expression for the correlation function (^(k, yu) with the shifted chemical 
potential /i as mentioned above 

G(k,/.)=L + ^A:^ + |^fc^) ' (40) 
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Now the conformational free energy of the chain under the influence of a force f can 
be calculated very easily by f ) = — lnG(0,/i;f) = — lnG(— i/^f, /i). Using the 
well-known thermodynamic relationship 

we express /i depending on its conjugate variable N and the force f. Note that the 
variable N is not the bare number of monomers A^^o since /i is a shifted chemical 
potential, but is proportional to Nq, indeed it is easily seen that N < Nq. Con- 
sidering only singular terms in k and neglecting terms of order this calculation 
yields: 

Substituting Eq. (H) into Eq. (g^) we get G{'k,N). According to Eqs. (|) and 
(H) (Rm) and (R5_) can be calculated from Eq. (|41|). Expanding in a power series 



for small forces to second order and again considering only most singular terms for 
small K, (Rm) becomes for Pf/n < 1 



(RS)[/] = + p'r4^ + ^43) 

The square root of the mean square elongation can be written according to Eq. (|43|) 



as 



. (Rj)[/]-(Rf)[0]~/ (44) 

which is a Hook-like law. 

Making the same approximation as mentioned above the root mean square end- 



to-end distance perpendicular to the force /, J (R^), decreases with / for < 1 
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which is contrary to a Gaussian chain 0. In particular: 

= - P'f-^^^ + Oif) (45) 

For / = the perpendicular end-to-end distance becomes exactly twice (Rm)[0]. 



6 Simulation 

In this section we briefly present the results from our Computer-simulations on a 
single polyelectrolyte chain. Again the aim is to get a force-size relationship for 
different values of the Debye-Hiickel screening length Ad- Therefore the monomers 
of the chain are located on the lattice points of a simple cubic lattice, i.e., the 
Kuhn segment length is equivalent to the distance between two neighboring lattice 
sites. Because of the fact that we only consider static properties of the chain, the 
algorithm of choice is clearly the pivot-algorithm, where one randomly chooses a 
link in the chain and then rotates this link together with the rest of the chain to a 
randomly chosen new orientation of the lattice PD|. Whether this configuration will 



be excepted or not is decided by a simple Metropolis-algorithm [^. 

We consider a chain of = 200 monomers with Debye-Hiickel interaction be- 
tween them. We made three different runs for the Debye-Hiickel screening length 
Ad = 5, 10 and 15 in units of the lattice constant. 

As the initial configuration we have chosen a totally stretched chain on the lattice. 
After 10^ pivot-steps, we defined the starting configuration for the further simu- 
lation. In the case without an applied force we made 1.6 * 10"^ pivot steps to get 
the end-to-end distance of the forceless reference chain. On the other hand in the 
case with an applied force we made 4 * 10^ pivot steps. To get sufficient statistical 
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Figure 2: Plot of the simulation data: Square root of the end-to-end distance parallel 
to the force minus the corresponding forceless case in double logarithmic scale 

Table 1 : Exponents for the force f in the force-size relationship for stretching of the 
chain parallel to the applied force 



Exponent 


Xd = 5 


Xd=10 


Xd = 15 


1. Regime 


0.67 


0.76 


0.84 


2. Regime 


0.4 


0.35 


0.34 



independent values for end-to-end distances we stored the configuration after every 
8000 steps. The results of the force-size relations are plotted in figure 2 and 3. 

On this double logarithmic scale we find for each screening length Xd two linear 
regimes. The change of the slope is clearly given in the region XoPf ~ 1, which is 
in excellent agreement with our theory. The slopes of the two regimes and therefore 
the exponent of the force / in the force-size relationship are given in table 1 for 
.y/(R|) and in table 2 for ^ (Ri). These values are calculated by a linear regression 
of the simulation data in figure 2 and 3. 
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Figure 3: Plot of the simulation data: Square root of the negative end-to-end dis- 
tance perpendicular to the force minus the corresponding forceless case in double 



logarithmic scale 



Table 2: Exponents for the force f in the force-size relationship for stretching of the 
chain perpendicular to the applied force 



Exponent 


Xd = 5 


\d = 10 


\d = 15 


1. Regime 


0.52 


0.56 


0.74 


2. Regime 


0.17 


0.14 


0.2 
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As one can see from table 1 and 2 the exponents of the first regime increase with 
increasing for both cases stretching parallel and perpendicular to the applied 
force f . According to our theoretical results we expect that the exponents for the 
first regime tend to one for larger Ad which is in good agreement with our data in 
table 1 and 2. For the second regime we find a drastically lower exponents in both 
tables. At this point it should be stressed that it is very difficult to get data for Xd 
greater than about 15, because in this case the apphed force f has to be so small that 
the transition from the first to the second regime lies within the numerical mistakes. 

7 Scaling considerations 

The analytical and numerical results suggest the following physical picture. At zero 
force, the polyelectrolyte chain is given by a chain of blobs. The blob size is entirely 
determined by the electrostatic properties (see Eq. (^31) ) . The low force regime 
(see fig. 4) can be viewed by the picture that the already greatly elongated chain of 
blobs becomes stretched. The chain size for zero force is determined by 



where we have ignored numerical prefactors. The fraction determined an effective 
step, or blob size, that is given by the charge z and the Debye screening length 
K,~^, i.e., oc z^/bXD■ FoT the latter equation we have assumed that the chain 
is weakly charged. Thus is consistent with the assumptions and the use of the 
random walk chain model. Remember that the variable is not the true degree of 
polymerization, but corresponds to the renormalized chain length via the relation 
(steepest descent Laplace inversion) A^ ~ l/(/^o + A^(0)). This chain of blobs can 




(46) 
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become elongated until the chain of blobs becomes fully stretched. In this low force 
regime the parts of the chain inside the blobs do not take part on the deformation 
process, (see fig. 4). This deformation process corresponds to the first regime in 
the simulation, i.e., for forces / < k^T/Xi) = fc- 

For larger forces (/ > fc) the blob size is no longer important, because the 
parts of the chain inside the blobs become deformed. This corresponds to a new 
Pincus regime, where the relevant blob size is now given by = keT//. The 
simulation clearly divides both regimes at forces Ad = k^T/f. Note that this fact 
has been already used in the analytical calculation above. The latter situation is 
very similar to the case considered by Rabin and Alexander, when the stretching 



of polymer brushes has been discussed |21|. Thus the blob size becomes diminished 



according to the idea pointed originally out by Pincus. This can be seen clearly 
from the simulations. The Pincus regime starts then at i? ~ Ad. 



8 Conclusions 

In the previous sections we analyzed the force-size relationship of a single polyelec- 
trolyte chain in a solvent. First we made theoretical considerations on this problem. 
Then we compared our results with computer-simulation data and some scaling con- 
siderations. In every three cases we find a transition from one to another regime 
at forces which are about XnPf ~ 1 where A^ = 1/k, the Debye-Hiickel screening 
length in the assumed Debye-Hiickel potential. The most important result of the 
present paper is the coupling between conformational degrees of freedom and the 
interactions in the elastic response of a single chain. The force is still proportional 
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f > f c (^\(S^(^6^6fio(b^ f > f c 

Figure 4: Blob-model of a polyelectrolyte chain with different external forces acting 
on the monomers (for details see text). 
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to the thermal energy fc^T, typical entropy elastic chains, but the single chain mod- 
ulus becomes strongly influenced by the interactions, via the Debye screening length 
Ad = K"^ 

Moreover the simulation data are in agreement with our theoretical predictions, 
that in the first regime the end-to-end distance of the chain depends linearly on the 
applied force in the long range limit of the Debye-Hiickel potential, which means 
in solvents with low salt concentration. In a subsequent paper we will extend this 



model and theoretical approach to the case of many crosslinked chains and make 
further predictions on the elasticity of polyelectrolyte networks. 
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